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The theoretical solution of the problem of the propagation of radio waxes through an 
ionized medium depends on Maxwell's well-known equations, together with the equation of 
motion of a free electron. This last equation, in general, contains a large number of terms. 
In addition to the purely electromagnetic terms, there are others expressing the interaction 
of the electron with the medium itself. The solution requires one to determine the complex 
index of refraction as a function of the applied radio frequency. The problem is even more 
involved, in some cases, because magnetohydrodynamical terms may enter into the ex- 
pression. Sound waves or shock waves in the medium, for example, can impress themselves 
on the complex index of refraction. As a result, dynamical and magnetohydrodynamical 
interactions can lead to the coupling of fields that otherwise would be independent. 



1. Introduction 

I wish, first of all, to extend a welcome on behalf 
of Dr. Ernest K. Smith and myself, cochairmen of 
this conference. My task this morning would have 
been much lighter if only two distinguished Soviet 
scientists, Prof. V. L. Ginsburg and Prof. A. V. 
Gurevich, had been able to accept our invitation. 
Their joint paper [1960] on ''Nonlinear Phenomena 
in a Plasma Located in an Alternating Electromag- 
netic Field" is a classic. It reviews and extends the 
basic principles and applies them to an examination 
of general nonlinear problems. 

Their paper is far too long and detailed for me to 
review here. Perhaps I should just commend it to 
your attention and then sib down. But such brevity 
is hardly in the tradition of conference chairmen. 
Besides, I should like to discuss some relatively 
unfamiliar nonlinear processes that may occur in 
certain types of ionized plasmas. 

Various contributors to this conference have asked 
me to explain the qualifying adjective, "nonlinear," 
which defines the type of processes we are discussing 
here today. Let me say, first of all, that "nonlinear" 
does not signify "irregular" or "nonuniform." The 
mathematician, the physicist, and the engineer each 
has his own definition of "nonlinear." At first sight 
these explanations may seem to be very different, 
but they are actually only alternate ways of ex- 
pressing the same basic facts. 

dE 



To the mathematician, quantities such as E, 



d 2 E 



dt' 



--Tjj) etc., are linear in that they depend only on the 

first power of E. Differential equations involving 
these quantities have solutions of the form: 
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And if two solutions are simultaneously present, the 
complete solution is simply the sum of the two, or 

E=Mt)+Mt). 

Mathematicians call this the principle of super- 
position. 

Mathematical quantities involving the square, or 
higher power, or cross product of the foregoing, such 

as E 2 , I -tt ) > E -r-> etc., are nonlinear. Mathemati- 
cally, we recognize the nonlinearity because the com- 
plete solution is not the sum of two elementary 
solutions. Cross-product terms exist. The two dis- 
turbances react on one another. For example, 

E*=M)+2Mmt)+Mt), 

and so on. 

Physically, E may represent an electric field that 
varies with the time. Any number of electromag- 
netic fields may exist simultaneously in any medium. 
As long as the disturbances are linear, we may repre- 
sent their total effect as the sum of the individual 
disturbances. But if the phenomenon is nonlinear, 
the total effect may depend on the cross product of 
the individual disturbances. 

The engineer has a practical definition of nonlinear 
phenomena. As long as a simple filter can separate 
the different superposed periodic disturbances passing 
through a circuit, they are linear. But when the 
two or more disturbances have interacted in such a 
way that they cannot be so separated, then the 
phenomenon is nonlinear. Under certain circum- 
stances where he wishes to induce artificial mixing, 
the engineer will use a nonlinear device such as a 
rectifier or detector. 

Many of the differential equations of mathematical 
physics are linear in the first-order. They may 
contain, however, certain nonlinear terms that are 
negligible unless the variable is exceptionally large. 
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Consider the following equation: 

rPE 
a ^+u 2 E+t3E 2 =0, 

where /3 is small. When /3 is zero, this equation 
has the solution: 

E=A sin <jit. 

But when we substitute this into the original 
equation, the term in E becomes 

(3E 2 =pA 2 sin 2 cot=^ (1-cos 2a>t). 

In other words, the perturbations have caused the 
harmonic of frequency 2co to appear. 

To obtain a general solution, we may employ 
the complex Fourier series: 



E= X) A n e iunt . 

n = —<=o 

Substituting this into the differential equation, we 
get 

Carry out the double summation by first summing 
over pairs of m and m f such that 

m -\-m' =n. 
Then we can write 

n L m J 

Hence we must have 



=0. 



A n a> 2 (n 2 -l)+pJlA m A n _ m =0. 



This infinite set of simultaneous equations is, of 
course, nonlinear in the A n 's. The solution contains 
higher harmonics, though they may have progres- 
sively smaller amplitudes. 

The most familiar nonlinear ionospheric process 
has received the name Luxembourg Effect, because 
it first showed up in connection with radio Luxem- 
bourg. The Luxembourg Effect is a particular 
example of a general phenomenon that might be 
more accurately termed "The Interaction of Acoustic 
and Electromagnetic Waves in an Ionized Medium." 
The medium may or may not contain a magnetic 
field. 

In the Luxembourg Effect, a radio transmitter of 
high power produces an intense electromagnetic 
field in the ionosphere. If an audio frequency 
modulates the carrier, the electron velocity field 
responds to this modulation. In effect, the tem- 
perature of the electron gas fluctuates by an ap- 
preciable amount, according to the impressed audio 
frequency. This temperature variation may also 



be described as a pressure variation, a true acoustic 
wave. I shall refer to the initial modulated wave 
and the acoustic disturbance it produces in the 
medium as the "unwanted wave. 7 ' 

The various propagation constants that depend 
on temperature or pressure of the electron gas in the 
medium vary with the frequency of the acoustic 
disturbance. If the density of the gas is high 
enough to make the medium dissipative, the imag- 
inary part of the complex index of refraction also 
varies. A second radio wave, termed the " wanted 
wave," traversing the medium, will suffer from this 
variable absorption. Hence the wanted wave ac- 
quires some of the modulation of the unwanted 
wave. Since this superposed modulation constitutes 
an interference with the wanted wave, it is undesired 
or " unwanted." Thus the terminology. 

The theory of this wave interaction was first 
given by Bailey and Martyn [1934], who pointed out 
the importance of magnetic fields, especially for 
radio waves near the gyrofrequency. The formulas 
they derived, however, were oversimplified, in that 
they linearized the equations describing the inter- 
action. As a result, their formulas do not properly 
display the phase lag of the acoustic disturbance 
relative to the modulation by an amount dependent 
on the audio frequency. Menzel and Layzer have 
developed a more detailed theory which will be 
presented at this conference. 

Experiments have shown that the Luxembourg 
variety of wave interaction rarely leads to cross 
modulation in excess of a few percent. The amount 
is very sensitive to a number of parameters. In the 
medium itself the significant variables are: the 
electron density and temperature, the frequency of 
dissipative electron collisions, and the magnetic 
field that fixes the electron gyrofrequency. The 
phenomenon further depends on the power and 
polarization, in addition to the radio and audio 
frequencies of the unwanted wave. The frequency 
and polarization of the wanted wave are also sig- 
nificant. The angle of incidence on the ionosphere 
determines the degree of penetration. Since the 
unwanted acoustic disturbance is localized, the 
greatest interaction will occur when the wanted 
waves travel through a significant portion of the 
disturbed volume. 

Although the original theory and observations 
refer to the interaction of two separate and inde- 
pendent waves, examination of the basic equation 
shows that a single modulated wave in the medium 
may interact with itself. As before, the fluctuating 
field causes an acoustic disturbance in the medium. 
This results in a variable attenuation coefficient, 
which reacts on the original wave. Since the phase 
lag is a function of the acoustic frequency, the 
degree of attenuation will similarly vary with the 
acoustic frequency. In effect, some frequencies will 
experience demodulation whereas others will be 
excessively modulated. 

The foregoing phenomenon, termed "self inter- 
action," indicates that any acoustic disturbance — 
of whatever origin — can impress its modulation on 



a wave traversing the volume. Rockets, missiles, 
explosions, collisions of solar ion clouds with the 
upper atmosphere, etc., represent potential sources 
of such acoustic waves. The question requiring 
study is the magnitude of the interaction, not 
whether such interaction occurs. More specifically, 
of the various kinds of interaction, are there any 
that could be studied experimentally? 

Certain difficulties are obviously present, resulting 
from fundamental differences in acoustic fields of 
radio origin on one hand and of mechanical origin on 
the other. For example, since the velocity of propa- 
gation of a radio wave greatly exceeds that of sound, 
the acoustic field induced by a modulated electro- 
magnetic wave will be in phase over a large volume 
of space. Hence the interaction with the wanted 
wave will likewise be in phase over a large volume of 
space. 

A mechanical disturbance, on the contrary, prop- 
agates with a much smaller velocity. Hence a 
wanted wave, progressing through the perturbed 
medium, will suffer high attenuation in certain zones 
and low attenuation in others. The effect resembles 
that of interference on optical frequencies. The 
problem is complex and the total effect on the wanted 
wave depends on the character and the extent of the 
basic disturbance. 

The details of the interaction require precise analy- 
sis. The limited volume of the mechanical disturb- 
ance may well be offset by the intensity of the dis- 
turbance itself. In a region where shock waves may 
exist, numerous nonlinear effects may occur, com- 
plicating the transport of electromagnetic waves 
through the medium. The complex (dissipative) 
index of refraction will probably contain significant 
velocity-dependent terms of nonlinear character. 
The effect on a radio wave transversing the medium 
may still be appreciable. Magnetohydrodynamic 
compression of the gas will significantly alter the 
magnetic field, with attendant complication of the 
mathematical problem. 

M. Cutolo has demonstrated experimentally the 
existence of a phenomenon sufficiently different from 
the Luxembourg Effect to deserve a special desig- 
nation. He refers to it as the "detection effect/ ' 
since it depends on nonlinear ionospheric properties 
for its occurrence. I prefer to call it the "Cutolo 
Effect/' after its discoverer. I shall describe it 
briefly, since Cutolo plans to devote his paper to 
more conventional nonlinear problems. 

Cutolo directs a beam of pulsed, modulated VHF 
waves upon the ionosphere. Of significance is his 
use of the gyro frequency as the modulation fre- 
quency. In the Luxembourg terminology, this 
pulsed, modulated, VHF wave, further modulated 
at audio frequencies to carry intelligence, represents 
the "unwanted wave." Cutolo reports that a sec- 
ond HF wave, the wanted wave, reflected from the 
ionospheric region traversed by the intense, focused 
VHF signal, will have the audio signal transferred to 
it. The amount of this cross modulation is very 
sensitive to the exactness of match between the gyro 
and modulating frequencies. Cutolo has found that 



the effect occurs for both vertical and oblique inci- 
dence of the wanted wave. 

I have not seen, up to this time, any quantitative 
physical explanation of the Cutolo effect. An un- 
wanted wave on the gyrofrequency would produce 
a large interaction of the Luxembourg type. But 
a VHF wave, modulated at the gyrofrequency, is 
not equivalent to a pure electromagnetic wave on the 
gyrofrequency. As is well known, such a wave can 
be expressed as an unmodulated carrier and (wo 
unmodulated sidebands separated from the carrier 
by the gyrofrequency. To account for the Cutolo 
effect we must suppose that some unknown iono- 
spheric nonlinearity so mixes the carrier and side- 
hands as to release the true gyrofrequency. It is 
difficult to see, though, how such interactions occur 
because of the very high frequencies employed. 
( )utolo used a carrier of the order of 50 Mc/s. 

The following analysis will serve as an elementary 
introduction to ionospheric theory for younger 
scientists. The paper also considers various non- 
linear effects. An electron of mass ra* and charge 
— 6,, moving with vector velocity v^ in an electro- 
magnetic field E and H and external acceleration F, 
obeys the well-known equation of motion: 



Mi 



dt 



,E— 2v,xH+wi,F. 



(1) 



This equation applies along with Maxwell's equations 
to define the physical state of the medium and its 
interaction with the field. 

The field F represents several varieties of mechani- 
cal actions. In a macroscopic problem involving 
the entire atmosphere, F will include the gravi- 
tational acceleration g. It also represents the 
collisional effects of nearby molecules, which may 
act as a sort of resistance proportional to the velocity 
and collision frequency v. As we shall note later on, 
this collisional term may be anisotropic, because of 
pressure gradients in the gas. 

Maxwell's equations are: 



1 TT TT * <^E . 47T T 

curl H=vxH -— H J, 

c ot c 



curl E=vxE= — 



c dt' 



div H=v • H=0, 



div E=V • E=4ttp c //c, 



(2) 



(3) 



(4) 



(5) 



wherein J is the electric current, p € the density of 
electric charge, k the dielectric constant and /x the 
magnetic permeability. Simultaneous solution of 



these equations leads to 



V 2 E=^ *J 



+VV-E, 



V 2 H= 



UK 



£> 2 H 

di 2 ' 



(6) 



(7) 



When p e is constant, the last term of (6) vanishes 
and we have the ordinary wave equations. For a 
plane wave, polarized in the zx plane and traveling* 
along the 2-axis, we have 



E=iE Q e ±u 



(8) 



where co is the circular frequency and c/n the veloc- 
ity of light in the medium. Thus n is the index of 
refraction. The quantity n may be complex. The 
unit vectors i, j, k are the customary Cartesian set. 
This equation satisfies (6) when 



n 2 =fxK. 



(9) 



For nonmagnetic media, ju=l. Hence k=1 for 
free space. Part of our problem involves the 
determination of the effective k when the space 
contains electric currents and charges. 
Return now to (1) in the form: 

d\ ■ e 

rrti -jt = — **E — - \xU—m i py l . (10) 

The first term on the right represents electric forces 
and the second, magnetic forces. The third term, 
somewhat schematically, represents the resistance 
that the charged particle encounters as it moves 
through the medium. The quantity v is the number 
of collisions per second experienced by the moving- 
ion. For convenience we shall drop the subscript i. 
The electrons, because of their great mobility, are 
the major contributors to the current density J. 

Thus: 



J = _AT €V , 



(11) 



where N is the number of electrons per cm 3 . Under 
different assumptions we shall calculate the varia- 
tions of v and J with the time. We shall then sub- 
stitute the J, so calculated, into Maxwell's equation 

(2). 

Case I. H=0, v=0, E=const:v ==— cEt/m. (12) 
Case II. E = 0, *>=0, H=const. (13) 

Here take the derivative of (10) with respect to the 
time. The resulting equation reduces as follows: 



dP\ 
dt = 



: -i(^ H Mi) 2(vxH)xH - (14) 



Let v have two components \ p parallel to and v s 
perpendicular (German: senkrecht) to the field. The 
triple vector product reduces to: 



(vxH)xH=£Pv s 



so that: 



and 



dt 2 



~-coi\ s ; 



d 2 \ p 
dt 2 



0, 



v s =v 0s e iw ^; v p =v 0p ; w L =eH/mc. 



(15) 
(16) 



The particle executes a spiral path with the Larmor 
frequency co L . The radius, a, and pitch, p, of the 
spiral are 

a = V 0s /<a L , p=V 0p /u L . 



Case III. E=0, H=0, v^const; v=v e- pt , 
corresponding to simple damping. 
Case IV. E=E Q e iat , H=0, *>=0, ^| 
\=(eE /ioom)e iut 



(17) 
(18) 

(19) 



Case V. At this point, let us consider a more general 
problem for which we assume an electromagnetic 
wave traveling parallel to the 2-axis. Let E x and E y 
be the vector electric fields and H x and H y the vari- 
able magnetic fields associated with the radiation. 
Let the magnetic field of the ionosphere be of magni- 
tude H, in the direction of the propagation. 

The general equation of motion, (10), breaks up 
into equations for the three components: 



dv x e j-, 



m 



eH 

mc 



v y —vv x 



dv y e -rp eH 

TT= E y -\ V x — Wy 

dt m mc 



dv z 
dt 



--—vV z 



(20) 



Multiply the second equation by i and add to the 
first, employing (16). The result is 

■^(v x +iv y ) = (io) L —v) (v x +iv y )—^(E x +iE y ). (21) 



Maxwell's field equations become: 

dEyl bH x dE x ldH y 
dz c dt dz 



c dt 
dH y _K()E x 4irNeV x 



dz c dt c 

dH x k &E y 4irNeVy 
dz c dt c 



(22) 



Multiply the second and fourth of these equations by 
) and combine as before, to give: 

^ (H x +iH p ) = -2 | {Ex +iE a )+^^(v x +iv v ). 



Now introduce the complex quantities: 

v x -\- iv v =Ve ±ia( - t ~ nz ' c) 

E x +iE v =Ee ±i " u - n!/c) 

H x +iH u =He ±i ^ t - nzlc) 



(23) 



(24) 



where V, E, and Hare complex constants. With the 
aid of (24), (21) and (23) become: 

[i(—a L ±o)+v]V=-eE/m (25) 

inE=H (26) 

icmH= - kwE± liriNeV. (27) 

These equations contain four complex unknowns: 
E, H, V, and n. Of these, n is the only one uniquely 
determined. H and V are, as one would expect, 
merely proportional to E. We may solve directly 
for n 2 : 

■ , AwNe 2 1 

mo) o:-\-co L -fip 

the complex index of refraction. Let 

n=n r Tik, 

where n T and n k are real. Then: 



n 2 r -k 2 =l- 



AirNe 



(*2-f(j0L 



2n r Jc= 



mco (co = Fw z ,) 2 +j' 2 



mo) (wTwl) 2 +j 



The propagation factor becomes: 



±iu(t—nz/c) == .p±iu(t — n r z/c) — ukz/c 



(28) 

(29) 

(30) 
(31) 

(32) 



density, J. Then, substituting this back into 
Maxwell's equation (2), we derived, in effect the 
dielectric constant k, which by (9) is equal to the 
square of the complex index of refraction. Equation 
(2) thus takes the successive forms: 



■ „ 1 dE . 4ttJ 
curl H= ^t+ 

c at c 



dE 4?rA^v *dE icoKE 



c dt 



c dt 



(33) 



since v, E, and their derivations all have the same 
time dependence, e iwl . 

The physical significance of k becomes more 
meaningful, perhaps, if we follow the steps indicated 
by (33). Here we may use the complex form of (2), 



. „ 1 dE . 4ttJ 
curl H=- ^r+- 

c ot c 



{ 



± 



iuE 4wNeV 



\ p±io)(t — nz/( 



= ±ioo+4irNe 2 
iu . 47rA r e 2 



{*! 



mc i(— oi L ±a))-*rv 



:] 



fi e ±iu>(t-nz/c) 



= _L l0)K fi e ±iu(t-nz/c) 
C 



(34) 



where we have used (25) to express V in terms of E. 
Now, solving for k we get 



--7h 2 =\- 



1 



(35) 



CO 0) =F 0) L T IV 

as before. 

Now note some of the assumptions implicit in 
these derivations of k. In the strictest sense we 
should have a different velocity, \ t , for each electron. 
We would then have evaluated the average v and J 
as follows: 



N\=-e SVp-€ fv/(tO<fo, 



(36) 



The last term of the exponent indicates an amplitude 
that decreases exponentially with the distance. 
The effect results from the collisions, which transfer 
a certain amount of energy to the medium. 

At this point, note the basic philosophy behind 
these elementary derivations of the propagation 
parameters. We have solved the general equation 
of motion to determine the velocity of the electron 
under the influence of external fields of force, 
electromagnetic, magnetic, and collisional. Our 
first objective was the calculation of the current 



where j(v) is the velocity distribution function for the 
electron gas. 

The solution of this problem is far more difficult 
than appears at first sight. The function f=f(v), 
will generally depend on time and position. The 
Boltzmann transport equation is basic for the 
determination of/. 



f r f t+v . vf+a . VrJ+s=0 . 



(37) 



In the above, a is the acceleration and v» the gradient 
in velocity space. The quantity, S, is the collision 
integral, which depends on the rate of change in the 
distribution function resulting from collisions. This 



term takes account of the creation or destruction of 
electrons in a given range of velocity. To calculate 
S we must know the physical details of both elastic 
and inelastic collisions. 

The first three terms of (37) are vector abbrevia- 
tions for the purely mathematical concept of the 
total derivative of a function, /, that depends on the 
time, on the coordinates, and the velocity components. 
In rectangular coordinates, we may write 



d 
dt 



5/ bfdx Zfdy dfdz 
bt^bx dt^by dt^~bz dt 



J(X, y, Z, U, V, co, 0=5t+$I ~/7+^7. 7/7 + 



^du ^dv Md. 
du dt dv dt dco dt 



wherein u, v, and to are the three velocity compo- 
nents. 

Let v' and v be the vector velocities of the electron 
prior to and after the collision. Let \[ and \ x be the 
initial and final velocities of the colliding particle. 
Let j(y')dy' and j(\)d\ symbolically represent the 
number of electrons within the ranges d\' and d\, and 
let F(\[)d\[ and Fiy^dVi similarly represent the 
number of colliding particles in ranges d\[ and d\\. 

The relative velocity of collision becomes: 



\— Vi and u=\\— vj 



(39) 



Let cr(u, 6) be the differential scattering cross section 
describing the probability that a colliding electron 
will be deflected through an angle 6, during collision. 
Define the solid angle rfO=sin dddd<t>. Then we may 
express S as follows: 

-ST — jjdvdQaiu, d)u[f(y)F(y 1 )-f(^)F(V 1 )}. (40) 

The problem of calculating S and / depends on a 
precise knowledge of the target cross section, a, as a 
function of u and 6. We shall not consider it further 
here, except to note that our arbitrary representation 
of the collisional effects by the term mv\ in (10) was 
naive, to say the least. It was particularly an over- 
simplification in that we assumed v to be independent 
of v, in order to obtain a linear equation. 

We should have been much more realistic if we had 
assumed v varied with the velocity. Certainly the 
collision frequency will depend on the temperature 
of the electron gas. This temperature, in turn, will 
vary with the impressed electric field. And if the 
field is modulated, the modulation frequency will 
cause v to vary at that frequency. Thus the absorp- 
tion coefficient, k, (31), should also vary with the 
time. It is this variable absorption that leads to 
cross modulation, the Luxembourg effect, gyro- 
interaction and other nonlinear responses, some of 
which we shall be discussing during this conference. 

But additional, still less obvious assumptions exist 
in our oversimplified theoretical discussion. Max- 
well's equations apply, of course. We have ignored 



possible reactions of the ionospheric plasma with the 
terrestrial magnetic field. I refer particularly to 
phenomena lying in the realm of magnetohydro- 
dyn amies. We can at least infer the character of 
some of these problems. 

Rewrite (1) in the general form: 



dVi ^ 



^xH 



•m j F, 



(41) 



wherein a positive e t refers to an ion and a negative 
€i to an electron. We can take averages over a unit 
volume as follows: 






(42) 



wherein p is the mass per unit volume. These sum- 
mations imply, when relevant, integration over the 
velocity-distribution functions previously discussed. 
The force-field F breaks up into two parts; external 
fields such as gravity and internal fields which result 
from colliding particles. One may write the latter 
force on a given particle, i, as follows: 



F int=S Fa, 



(43) 



where the summation is taken over all neighboring 
particles j. This term cancels out only when the 
medium is uniform. If electric currents are present, 
the term gives rise to forces simulated in our earlier 
treatment by the collisional term, m^v. This force 
can also be expressed, for the entire medium, as a 
sort of electrical resistance, R, in the form i?J. It 
also produces such major forces as pressure gradients. 
Tangential forces, in the presence of sheer, lead to 
the phenomenon of viscosity. In the presence of 
temperature gradients it can even cause thermo- 
electric effect. For example we may write 



E=B(J-9vT), 



(44) 



where 6 is a thermoelectric coefficient. In the strict- 
est sense both R and 6 will be tensors rather than 
scalars. 

As we sum (41) over the individual particles and 
employ Maxwell's equations we obtain the general 
equation of magnetohydrodynamic motion: 

P§ t = -Vp+pF exp +^-(vXH)XH+^(vXE) 

at 47T 47T 

<E+iE Vl .E- 4 -LA( EX H). (45) 



When II and E are constant or zero, the equatio 
reduces to that of ordinary hydrodynamics, less the 
second-order terms related to viscosity. The terms 
represent the forces arising from the electromagnetic 
fields. EXH is the Poynting vector, which repre- 
sents energy flowing out of the volume element. 
Hence the partial time derivative of EXH represents 
the forces caused by radiation pressure. 

Now multiply (31) by e/m and sum again. The 
result is 

^^=E+i(vXH)-£(J-0vT) (46) 



Now take the partial time derivative of (49) and 
derive a similar equation for E. The result is: 

2>E . 1 d , w „, c 2 BV k c) 2 E "I ,„. 

ar+- C ^ (vXH) -^L V2E ^^-W-E} (51) 

These equations are also subject to the equation of 
continuity, 

dp dp . 

_ = _ +v . vp= _ pV . y , 



or 



f?+V.(pv)=0. 



(52) 



where m is the mass of an electron and n the number 
of electrons per unit volume. In taking this sum 
we have allowed for the fact that electrons, by virtue 
of their small mass, are the main contributors to the 
electric currents. In passing we may note that this 
equation has several obvious limiting forms. It can 
reduce to Ohms law 



E=RJ, 



(47) 



or to the electromotive force induced in a conductor 
moving through a magnetic field 



E= 



-(vXH)/c. 



(48) 



Even (46) has its limitations. Acoustic waves 
or hydrodynamic waves can directly produce changes 
in J, by simple compression or expansion of the 
current-bearing regions. Moreover, in a conducting 
medium, H itself may change under magnetohydro- 
dynamic forces. Such changes in H produce fluc- 
tuations in J. As we have just noted, variations in 
J directly affect the dielectric constant. We must, 
therefore, set up the equations governing the varia- 
tions, not only of H but also of E. 

The following analysis is by no means complete. 
However, consider the steady-state situation, derived 
from (46) : 



E+(vXH)/c-#J=0. 



(49) 



Take the curl of this equation and reduce by means 
of the Maxwell relationships. We get : 



r/H 

~dt~~ 






+ v vH 



XT _ Uw ,C 2 £[~2T7 « d 2 H~| 



(50) 



For regions where the conductivity is high, /?— >0 
and we recover the ordinary equations of magneto- 
hydrodynamics. When the conductivity is low, 
divide through by 11. The electromagnetic equa- 
tions then reduce to the simple wave equations for 
E and H. In other words, for media of low conduc- 
tivity, the hydrodynamic and electromagnetic fields 
are independent. Otherwise the two fields are 
coupled. This means that, in effect, we can combine 
the purely dynamic terms into a nonlinear wave 
equation, in which the dielectric constant is a function 
of the dynamical variations. 

I shall not attempt to carry the analysis further. 
My primary objective has been to show that non- 
linear effects abound in this general problem. I 
should point out one additional fact. The electrical 
resistance R is closely related to the damping factor k. 
Since k, in turn, depends on both electrical and 
dynamical factors, the factor R introduces addi- 
tional nonlinearities. Detailed solutions depend on 
assumptions concerning the ba ic physical conditions 
in the medium. I have especially desired to empha- 
size the importance of dynamical or magnetohydro- 
dynamical effects on radio propagation. 
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